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Experiments on (Ga,Mn)As in the low-doping insulating phase have shown evidence for the pres- 
ence of an impurity band at 110 meV above the valence band. The motivation of this paper is to 
investigate the role of the impurity band in determining the magnetic correlations in the low-doping 
regime of the dilute magnetic semiconductors. For this purpose, we present results on the Haldane- 
Anderson model of transition- metal impurities in a semiconductor host, which were obtained by 
using the Hirsch-Fye Quantum Monte Carlo (QMC) algorithm. In particular, we present results on 
the impurity-impurity and impurity-host magnetic correlations in two and three-dimensional semi- 
conductors with quadratic band dispersions. In addition, we use the tight-binding approximation 
with experimentally-determined parameters to obtain the host band structure and the impurity- 
host hybridization for Mn impurities in GaAs. When the chemical potential is located between the 
top of the valence band and the impurity bound state (IBS), the impurities exhibit ferromagnetic 
(FM) correlations with the longest range. We show that these FM correlations are generated by the 
antiferromagnetic coupling of the host electronic spins to the impurity magnetic moment. Finally, 
we obtain an IBS energy of 100 meV, which is consistent with the experimental value of 110 meV, 
by combining the QMC technique with the tight-binding approach for a Mn impurity in GaAs. 

PACS numbers: 75.50.Pp, 75.30.Hx, 75.40.Mg, 71.55.-i 



I. INTRODUCTION 

The discovery of dilute magnetic semiconductor 
(DMS) materials is important because of possible spin- 
tronics device applications [3, 0, H, 0]- The electronic 
state of the alloy (Ga,Mn)As, where Mn substitutes 
for Ga, has been investigated by various experimental 
methods including transport measurements, optical and 
photoemission spectroscopy, and scanning-tunnelling mi- 
croscopy d H 0, B S 03 • In the low-doping regime 
(<C 1%), this prototypical DMS material is insulating 
with clear experimental evidence for the presence of an 
Mn-induced impurity band located 110 meV above the 
valence band. For more than 2% Mn doping, the dc 
transport measurements indicate that the impurity band 
emerges with the valence band leading to a metallic state. 
The redshift of the mid-infrared peak observed in opti- 
cal absorption in metallic samples has been attributed to 
inter-valence band transition instead of being due to an 
impurity band [1, [H|. Hence, the analysis of the exper- 
imental data suggest that (Ga,Mn)As exhibits two dif- 
ferent regimes at low and high impurity concentrations 
separated by an insulator- metal transition fllT ] . 

In this paper, our purpose is to develop a microscopic 
understanding of the low-doping insulating phase of the 
DMS material (Ga,Mn)As and to shed light on the role 
of the impurity band in producing the magnetic correla- 
tions. For this purpose, we present numerical results on 
the single- and two-impurity Haldane- Anderson model of 
transition-metal impurities in semiconductors (l2| using 
the Hirsch-Fye Quantum Monte Carlo (QMC) method 
[lU . We study the influence of the electronic state of the 
semiconductor host on the impurity bound state (IBS) 



and, in turn, on the magnetic correlations which develop 
between the impurities and between the impurity and the 
host electrons. We find that, when the chemical potential 
\i is located between the valence band and the IBS, the 
ferromagnetic (FM) correlations between the impurities 
exhibit the longest range. We also show that the inter- 
impurity FM correlations are generated by the antiferro- 
magnetic (AFM) correlations between the impurities and 
the host electrons. In addition, we use the tight-binding 
approximation to determine the host band structure and 
the impurity-host hybridization for the three ti g orbitals 
of Mn in GaAs. Using tight-binding parameters con- 
sistent with photoemission measurements, we obtain an 
IBS energy which is in agreement with the experimental 
value of 110 meV. Because of these results, we think that 
a microscopic understanding of the low-doping insulat- 
ing phase of (Ga,Mn)As is possible within the Haldanc- 
Anderson model. 

The nature of the magnetic correlations in the 
Haldane-Anderson model of transition-metal impurities 
in semiconductors was studied using the Hartree-Fock 
(HF) approximation [13 EH- It was shown that, when 
the chemical potential is located between the top of the 
valence band and the IBS located at loibs, long-range 
FM correlations develop between the impurities medi- 
ated by the AFM coupling of the valence electrons to 
the impurity magnetic moments. On the other hand, 
when the IBS becomes occupied, the spin polarization 
of the host split-off state cancels the polarization of the 
valence band, reducing the range of the FM correlations. 
The QMC calculations for the Haldane-Anderson model 
support the HF picture for the role of the IBS in produc- 
ing the FM correlations [lj|. The QMC results on the 
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range of the FM correlations are in agreement with the 
HF predictions, however, the HF approximation under- 
estimated the value of luibs- These results on the mag- 
netic correlations for < /i < loibs are different than 
in the metallic case, /i < 0, where the inter-impurity 
magnetic correlations exhibit Rudcrman-Kittel-Kasuya- 
Yosida (RKKY) oscillations with period determined by 
the Fermi wavevector kp- We note that the Haldane- 
Anderson model was also studied within HF by Krstajic 
et al. for DMS materials and by Yamauchi et al. [IH 
for hemoprotein. The role of IBS for DMS materials was 
also discussed by Inoue et al. [19| within HF. 

In this paper, we are particularly interested in how the 
magnetic properties are influenced by the host electronic 
state, which we model using two different approaches. 
First, we consider the simple case where the magnetic 
impurity has one orbital only and the host band struc- 
ture consists of quadratic valence and conduction bands. 
Here, in addition, we treat the impurity-host hybridiza- 
tion as a freely adjustable parameter. This single-orbital 
model is described in Section 2. A, and the QMC results 
for this case are presented in Sections 3 and 4 for two 
and three-dimensional semiconductors, respectively. We 
study the inter-impurity and impurity-host magnetic cor- 
relations, the induced electron-density around the impu- 
rity, and show how these quantities depend on the param- 
eters of the Haldane- Anderson model and the dimension- 
ality of the semiconductor. These results show that the 
IBS and the magnetic correlations depend sensitively on 
the model parameters. 

In order to develop a more realistic model for 
(Ga,Mn)As, we next use the tight-binding approxima- 
tion to obtain the band structure of the bare host GaAs 
and the hybridization matrix elements with the three ti g 
orbitals of Mn substituted in place of Ga. In this ap- 
proach, the tight-binding parameters required for calcu- 
lating the hybridization are taken from an analysis of 
the photoemission data on Mn 2p core level. This model 
is introduced in Section 2.B and the QMC results ob- 
tained for this case are presented in Section 5. These 
QMC calculations keep all three of the Mn t^g orbitals, 
and hence the multi-orbital effects are included except for 
those due to Coulomb repulsion between different Mn or- 
bitals. Within this approach, we obtain an ljibs which 
is close to the experimental value of 110 meV. We also 
show that the FM correlations between the impurities 
weaken as the IBS becomes occupied. These results em- 
phasize the importance of the IBS in the low-density limit 
of (Ga,Mn)As. 

Finally, we note that the approaches taken in this pa- 
per can be extended to the case of finite concentration 
of Mn impurities in order to study the insulator-metal 
transition and the metallic state of (Ga,Mn)As observed 
at higher Mn concentrations. 



II. IMPURITY MODEL 

The general model for describing transition-metal im- 
purities in a semiconductor host is given by Haldane- 
Anderson Hamiltonian [l2j , 

H = Y ( £ka ~ ' i ) C kaa C k« CT + Y ( Ed Z ~ ^) d ha d ^ 

+ Y (^ka^cL^+H-C-) 

k,o:,2,£,<7 

hi 

where c ] [ CQCT (ckao-) creates (annihilates) a host electron 
with wavevector k and spin a in the valence or conduc- 
tion bands denoted by a, and d\^ a is the creation 
(annihilation) operator for a localized electron at impu- 
rity orbital £ located at site i. The first term in Eq. (1) 
represents the kinetic energy of the host electrons, and 
the second term is the bare energy of the localized impu- 
rity orbitals, while the third term is due to the impurity- 
host hybridization. The last term represents the onsite 
Coulomb repulsion at the impurity orbitals. We note that 
here we are neglecting the Coulomb repulsion among the 
different impurity orbitals. The effects of the Hund cou- 
plings will be considered in a separate paper. As usual in 
Eq. (1), U is the onsite Coulomb repulsion, fi the chem- 
ical potential and is the bare energy of the impurity 
orbital £. In addition, the hybridization matrix clement 
is 

l k„,: \u„,< k " ■ (2) 

where is the coordinate of the impurity site j and 
Vka.f is the value when the impurity is located at the 
origin. We use the Hirsch-Fye QMC technique to study 
the Haldane- Anderson impurity Hamiltonian, Eq. (1), 
for the the single- and two- impurity cases [1 31 ] . 

In this paper, we perform the QMC calculations us- 
ing two different approaches for incorporating the host 
band structure and the impurity-host hybridization. In 
Sections 3 and 4, we consider a simple model where the 
impurity site contains a single impurity orbital and the 
semiconductor bands have quadratic dispersion, while 
in Section 5 we use the tight-binding approximation for 
treating the three ti g orbitals of Mn impurities in GaAs. 
In the remainder of this section, we describe these two 
different approaches for modelling the transition-metal 
impurities in semiconductors. 

A. Single-orbital case 

In Sections 3 and 4, we consider the simple case 
where the impurity site contains a single orbital in two 
and three-dimensional semiconductor hosts, respectively. 
Hence, in this case Eq. (1) reduces to Eq. (1) of Ref. [lj|. 
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FIG. 1: (Color online) Schematic drawing of the semicon- 
ductor host bands e£ (solid curves) and the impurity bound 
states (thick arrows) obtained with HF in the semiconductor 
gap. The dashed line denotes the chemical potential fi. 

Furthermore, in these sections, we assume that the host 
band structure consists of one valence (a = v) and one 
conduction (a — c) bands with quadratic dispersions 
given by 

ek,„ =-D(k/k ) 2 (3) 
e k , c = D(k/k ) 2 + A G . (4) 

Here, D is the bandwidth, ko the maximum wavevector 
and Ag the semiconductor gap. The energy scale is de- 
termined by setting D = 12.0. Figure 1 shows a sketch 
of the host band structure. For the Coulomb repulsion 
we use U — 4.0, and set the bare value of the impurity 
orbital energy to Ed = fi—U /2, so that the impurity sites 
develop large magnetic moments. We report results for 
semiconductor gap A G = 2.0, and inverse temperature 
(3 = l/T from 4 to 32. We also use a constant V for 
Vka,£j an d treat V as a free parameter. 

For the single orbital cases of Sections 3 and 4, we show 
QMC results on the equal-time impurity-impurity mag- 
netic correlation function (Mf Mf ) for the two-impurity 
Haldane-Anderson model. Here, the impurity magneti- 
zation operator is given by 

Mf = 4 r dn - 4i<kl, ( 5 ) 

and the fermion creation and annihilation operators act 
at a single orbital at the impurity site. In addition, 
we present results on the impurity-host correlation func- 
tion (M z m z (r)) for the single-orbital and single-impurity 
Haldane-Anderson model. Here, the host magnetization 
at a distance r away from the impurity site is given by 

m z (r) = J2 (4 T W c «tW ~<li(r)c a l(r))- (6) 

a—v,c 

For the metallic case, the correlation functions (Mf Mf ) 
and (M z m z (r)) were previously studied by using QMC 
[HI, 0i HH, H2|. In the single impurity case, we also 



present QMC data on the square of the impurity moment, 
((M z ) 2 ), and the impurity susceptibility x defined by 

X = / dr{M z {r)M z {Q)). (7) 
Jo 

The effects of IBS are clearly visible in these quantities. 
We also present QMC data on the charge distribution 
of the host material around the impurity. These QMC 
results are obtained using Matsubara time step At = 
0.25. At j3 = 16, (Mf Mf ) varies by a few percent as At 
decreases from 0.25 to 0.125. 



B. Tight-binding model for Mn in GaAs 

In order to construct a more realistic model of 
(Ga,Mn)As in the dilute limit, we use the tight-binding 
approach to calculate £ka and Vk a ,i£ of Eq. (1) for the 
Mn t 2g orbitals. In Section 5, we will present QMC data 
obtained using these results on e^ a and Vka,i£- 

In this approach, the tight-binding band structure £ka 
of GaAs host is obtained by keeping one s orbital and 
three p orbitals at each Ga and As site of GaAs with 
the zincblende crystal structure. We consider only the 
nearest-neighbor hoppings between the Ga and As sites. 
Figure 2 shows the resulting band energies Sk a versus k 
along various directions in the face-centered-cubic (FCC) 
Brillouin zone for pure GaAs. The Slater-Koster parame- 
ters [23[ which have been used for obtaining these results 
were taken from Ref. (24|. In Fig. 2, the top of the va- 
lence band is located at the T point, where the energy 
gap is 1.6 eV, which is consistent with the experimental 
value. This simple approximation, with nearest-neighbor 
hopping among the eight sp 3 orbitals, reconstructs the 
four valence bands reasonably. They are important for 
the impurity-host magnetic coupling in (Ga,Mn)As, be- 
cause the chemical potential is located near the top of 
the valence band. 

When a Mn impurity is substituted in place of Ga in 
GaAs, the five 3d orbitals of the Mn ion are split by the 
tetrahedral crystal field into the three-fold degenerate t 2g 
orbitals and the two-fold degenerate e g orbitals. Since 
the e g orbitals have bare energies which are much lower 
than the t 2g orbitals, we neglect the e g orbitals and keep 
only the three t 2g orbitals of Mn in Eq. (1). In order to 
calculate Vk a ,i? for the Mn t 2g orbitals, we assume that 
the sp 3 orbitals of Mn are the same as those of Ga, and 
this way the periodic boundary conditions are also sat- 
isfied. In other words, at the Mn impurity site, we only 
add the three Mn t 2g orbitals to the GaAs host. This is a 
simple model for a transition-metal impurity substituted 
into GaAs. However, this approach allows us to take into 
account the effects of the host band structure beyond the 
quadratic dispersion. Furthermore, in Sections 3 and 4, 
we treat the hybridization matrix element V as a free 
parameter. However, here, we perform the calculations 
using realistic parameters for the hybridization of the 3d 
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FIG. 2: (Color online) Band structure of GaAs within the 
tight-binding approximation. 

orbitals with the host semiconductor bands. We also note 
that in Ref. [25[ the tight-binding and QMC techniques 
have been combined to study how the crystal structure 
of ZnO host affects the magnetic properties when Mn im- 
purities are substituted. However, in these calculations 
only one of the Mn 3d orbitals is taken into account and 
the multi-orbital effects are not included. 

Within the tight-binding approximation, the hy- 
bridization matrix element Vk Q .£, Eq. (2), is obtained 
from 

*W = (tfka|floN>, (8) 

where \ipua) is the Bloch eigenstate of the host electrons, 
\ipt-) the localized orbital eigenstate of an impurity lo- 
cated at the origin, and Hq is the Hamiltonian of the 
system with the Coulomb repulsion turned off. In order 
to evaluate Vka,£, it is necessary to determine the values 
of the Slater-Koster parameters between the sp 3 orbitals 
of GaAs and the d orbitals of Mn, which are denoted by 
the notation (spa), (pda) and (pdn) 23]. In the follow- 
ing, we determine the value of (pdir) by the general equa- 
tion (pdir) = (pda)/ (-2.16) \2(y\. We also set (sda) = 
for the following reason: When the chemical potential 
is near the top of the valence band, the hybridization 
matrix elements near the T point become important for 
the magnetic couplings. However, the s orbitals do not 
contribute to the three degenerate valence bands at the 
top of the valence band [27| ■ This is why we ignore the 
contribution of the s orbitals to hybridization and set 
(sda) = 0. The experimental estimate for the remain- 
ing Slater-Koster parameter is (pda) = —1.1 eV, which 
is obtained from the photoemission data on the Mn 2p 
core level with a configuration-interaction analysis based 




FIG. 3: (Color online) Hybridization matrix element V\t, a ,xy 
of the £ = xy orbital of a Mn impurity with (a) the valence 
and (b) the conduction bands of GaAs versus k along various 
directions in the FCC Brillouin zone obtained using the tight- 
binding approximation. 

on a cluster model [28]. In Section 5, we will present 
QMC data obtained using (pda) = — 1.14 eV yielding an 
ujibs — 100 meV, which is in reasonable agreement with 
the experimental value of 110 meV for (Ga,Mn)As in the 
low-density limit. 

Figure 3 displays tight-binding results on Vk Q ^ versus 
k for the £ = xy orbital, and (a) the valence and (b) 
conduction bands along various cuts in the FCC Bril- 
louin zone. The values of Vk a ,£ at the top of the valence 
and at the bottom of the conduction band are particu- 
larly important for the IBS. We note that Vk a ,f=xy takes 
large values near the T point for the top valence bands, 
while it is weaker for the conduction bands. In fact, the 
hybridization of the xy orbital with the lowest-lying con- 
duction band vanishes at the T point. We also notice that 
Vk Q ^ can be discontinuous around the high-symmetry 
points. In Section 5, we will present QMC data obtained 
using these values for e^ a and Vka.{- In addition, for 
the bare energies of the impurity orbitals, we will use 
Ed$ = [i—U/2 with U = 4.0, so that large magnetic mo- 
ments develop at the impurity sites. The QMC results 
do not depend sensitively on small changes in the values 
of Ed£ and U. 

After obtaining £ka and Vkc*,^ within the tight-binding 
model described above, we perform the QMC calculations 
keeping all three Mn t2 g orbitals. Hence, the QMC results 
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which will be presented in Section 5 include the multi- 
orbital effects except for the Hund coupling. However, 
we will present these QMC results only for the £ = xy 
orbital in order to make comparisons with the the single- 
orbital cases of Sections 3 and 4. In particular, defining 
the magnetization operator for £ = xy as 



M?=d\^d ia -d\ a d ia , 



(9) 



we will present data on ({M z ) 2 ) and the magnetic sus- 
ceptibility x, Eq. (7), for the xy orbital within the single- 
impurity Haldane- Anderson model. For the two-impurity 
case, we will present data on the inter-impurity magnetic 
correlation function (Mf Mf ) between the xy orbitals lo- 
cated at the impurity sites Ri and R2. In addition, we 
will discuss the magnetic correlation function between 
the magnetic moment of the xy orbital and the host elec- 
trons (M 2 m z (r)), where the host magnetization operator 
at lattice site r is 



rat 



C raJ, C raJ.) 



(10) 



with a summing over the eight semiconductor bands, for 
the single impurity case. In Section 5, we will see that the 
experimental value of loibs is reproduced reasonably by 
combining the tight-binding approach with the QMC cal- 
culations. In addition, we will see that the magnetic cor- 
relations weaken as the IBS becomes occupied, in agree- 
ment with the results of Sections 3 and 4. 

In the following sections, we will find that the quanti- 
tative results depend on the dimensionality and the band 
structure of the host material. In Sections 3 and 4, we 
show results for the two and three dimensional host ma- 
terials with simple quadratic band dispersions. Then, in 
Section 5, we discuss the case for a GaAs host using the 
tight binding approximation. 



III. QMC RESULTS FOR A 2D 
SEMICONDUCTOR HOST 
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FIG. 4: (Color online) Inter-impurity magnetic correlation 
function (MfMI) vs koR at various j3 for (a) /1 = —6.0 and 
(b) \i = 0.1 for the two- impurity Haldane- Anderson model. 
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In this section, we show results for the 2D case. Here, 
the density of states of the pure host, po, is a constant 
with a sharp cutoff at the semiconductor gap edge, which 
leads to stronger impurity-host coupling compared to the 
3D case. Here, we present results for hybridization A = 
irp V 2 varying from 1 to 4. 



A. Magnetic correlations between the impurities 

Figures 4(a) and (b) show (MfMf ) versus koR, where 
R is the impurity separation, for the two-impurity An- 
derson model for half-filled metallic (p = —6.0) and semi- 
conductor (p = 0.1) cases. In the metallic case, (M-fMf) 
exhibits RKKY-type oscillations with both FM and AFM 
correlations depending on the value of k$R. On the other 



FIG. 5: (Color online) Impurity magnetic-moment square 
{(M z ) 2 ) vs p at various (3 for A = 1, 2 and 4 for the single- 
impurity Haldane- Anderson model. 



hand, for p — 0.1, we observe FM correlations of which 
range increases with (3. 

The remainder of the data shown for the 2D case in 
this section are for the single-impurity Anderson model. 
In Fig. 5, we show results on ((M z ) 2 ) vs p for various 
values of A. As T is lowered, a step discontinuity de- 
velops in ((M z ) 2 ) near the gap edge. The location of 
the discontinuity coincides with the location of IBS de- 
duced from data on (M-f M|), the impurity single-particle 
spectral weight A(u>), and the inter- impurity susceptibil- 
ity X12 discussed previously [16| . Hence, Fig. 5 shows 
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FIG. 6: (Color online) T\ vs fi at various j3 for A = 1 for 
the single-impurity Haldane- Anderson model. 



that the local moment increases rapidly as the IBS be- 
comes occupied. Here, we also observe that the moment 
size decreases with increasing hybridization, as expected. 
Figure 6 shows T\ vs /j, for A = 1, where we observe that 
a step discontinuity develops at the same location as in 
((M 2 ) 2 ) shown in Fig. 5(a). For /i ^ ujjbs, T\ decreases 
with decreasing T due to the screening of the impurity 
moment by the valence electrons. However, for /i ^ ujibs, 
the impurity susceptibility exhibits free-moment behav- 
ior in agreement with the role of the IBS discussed above. 
We note that determining the location of the IBS from 
A(lu) is costly in terms of computation time. For this 
reason, in the single-impurity case, it is more convenient 
to determine uoibs from data on ((M z ) 2 ) versus \i. In 
the remaining sections, we will use this approach for de- 
termining UJIBS- 

Within the HF approximation and in 2D, the range of 
the FM correlations is given by (IGttpoujibs)^ 1 ^ 2 for < 
(i < ujibs- This implies that the range decreases with 
increasing ojibs and, hence, A. We find that the QMC 
results on the maximum range of the FM correlations 
are in quantitative agreement with the values from the 
HF calculations [l6[. However, we also find that the HF 
approximation underestimates the value of ujibs- F° r 
example, for A = 1, HF yields ujibs — 0.017, which is 
about an order of magnitude smaller. Hence, within HF, 
the long-range FM correlations are restricted to a narrow 
range of fi. 



B. Impurity-host correlations 

In this section, we discuss the magnetic correlations 
between the impurity and the host. In addition, we show 
results on the induced charge oscillations around the im- 
purity site. 

Figures 7(a)-(c) show the impurity-host magnetic cor- 
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FIG. 7: (Color online) Impurity-host magnetic correlation 
function s(r) vs kor at various (3 for (a) \i — —6.0, (b) /i = 0.1, 
and (c) /i = 0.5 for the single-impurity Haldane- Anderson 
model. 



relation function s(r) defined by 



s(r) 



li^k^r 



no 



(M z m z {r)) 



(11) 



versus kor for the single-impurity Anderson model for 
/i = —6.0 (half-filled metallic), /i = 0.1 (semiconductor 
with IBS unoccupied), and fi = 0.5 (semiconductor with 
IBS occupied). Here, hq is the electron density and r is 
the distance from the impurity site. For fi — 0.1, the 
coupling between the impurity and host spins is AFM 
for all values of k^r, while, for [i = —6.0, s(r) exhibits 
RKKY-type 2kp oscillations. We also note that, for the 
metallic case, the magnetic correlations saturate before 
reaching f3 — 32. Comparison of Fig. 7(b) and Fig. 4(b) 
for \i — 0.1 show that the AFM impurity- host coupling 
produces the FM correlations between the impurities. 

As seen in Fig. 7(c), when (j, is increased to 0.5, the 
AFM correlations between the impurity and host spins 
become much weaker. This is because the IBS becomes 
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FIG. 8: (Color online) S(kor = 25) vs /i at various j3 for the 
single-impurity Haldane- Anderson model. 



occupied for \i > 0.1. Within the HF approximation 
0, Gil; the FM interaction between the impurities is 
mediated by the impurity-induced polarization of the va- 
lence electron spins, which exhibit an AFM coupling to 
the impurity moments. The impurity-host hybridization 
also induces host split-off states at the same energy as 
the IBS. When the split-off state becomes occupied, the 
spin polarizations of the valence band and the split-off 
state cancel. This causes the long-range FM correlations 
between the impurities to vanish. These QMC and HF 
results emphasize the role of the IBS in determining the 
range of the magnetic correlations for a semiconductor 
host. 

The total magnetic coupling of the impurity magnetic 
moment to the host is obtained from 



S(k r) 



k r 



d(k r')s(r') 



(12) 



Figure 8 shows S(kor = 25) vs fj, for A = 1 at various 
f3. Here, we see that the impurity becomes magnetically 
decoupled from the host when h^ujibs- 

Next, in Fig. 9, we show the modulation of the charge 
density around the impurity. Here, we plot p(r) vs kor, 
where p(r) is defined by 



p{r) 



^ 2 ° (n a (r) - n a (oo)) 



no 



(13) 



with n a (r) = J2<j( c acr( r ) c av (»")). For the metallic case 
of /i = —6.0, we observe long-range RKKY-type of os- 
cillations in p(r). When /i = 0.1, the charge density 
around the impurity is depleted up to koR ~ 20 at these 
temperatures. This depletion represents the extended 
valence hole which forms around the impurity. However, 
for /i = 0.5, the induced charge density decreases signifi- 
cantly as T is lowered, because now the IBS is occupied. 
We next integrate p(r), 



P(k r) 



k r 



d(k r')p(r'), 



(14) 
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FIG. 9: (Color online) p(r) vs kor for A = 1 and (a) /i = 
—6.0, (b) 0.1 and (c) 0.5 for the single- impurity Haldane- 
Anderson model. 



and plot P(kor — 25) vs [i in Fig. 10. We observe that 
the total charge density around the impurity is most de- 
pleted when /i^Sljibs, which is due to the valence 
hole induced around the impurity. In the metallic case, 
the induced charge density is oscillatory and has a long 
range as we have seen in Fig. 9(a). As /i approaches 
the gap edge, the electron density around the impurity 
is depleted. However, when fi^tujiBS: w e see that this 
depletion is cancelled by the extended charge density of 
the split-off state. 



IV. QMC RESULTS FOR A 3D 
SEMICONDUCTOR HOST 

In this section, we discuss the three-dimensional case, 
where the hybridization parameter irV 2 N(0) vanishes at 
the gap edge because of the vanishing of N(0) of the 
pure host. Hence, the impurity-host coupling near the 
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FIG. 10: (Color online) P(kor — 25) vs p at various (3 for 
the single-impurity Haldane- Anderson model. 



gap edge is much weaker compared to the 2D case. Con- 
sequently, the FM correlations between the impurities 
have a shorter range. We find that the dimensionality of 
the host material strongly influences the magnetic corre- 
lations. In particular, we see that the IBS does not exist 
in 3D unless the hybridization is sufficiently large. 

Here, we define hybridization as A = irp^V 2 where p$ 
is the density of states when the valence band is half- 
filled, Pq = kl / (A^/2ii 2 D) . This choice allows us to use 
comparable values for the hybridization matrix element 
V when we compare the 2D and 3D results. In the fol- 
lowing, we present results for A = 1, 2 and 4. 

Figures ll(a)-(c) show (Mf Mf ) vs k$R at various val- 
ues of p for A = 1, 2 and 4. These results are for 
p = —6.0, p = —1.0 and p = 0.0. We observe that, 
for A = 1, the FM correlations between the impurities 
weaken as p approaches the top of the valence band. On 
the other hand, for A = 2 and 4, the FM correlations are 
stronger for p = 0.0. This is because, for A = 1, the IBS 
does not exist in a 3D host, as we will see later in Fig. 
13, which shows results on ((M z ) 2 ) versus p. 

Figure 12 shows the impurity-host magnetic correla- 
tion function s{r) vs kor for the single-impurity case. In 
3D, s(r) is defined by 



s(r) 



47r(fcor) 
n 



■(M z m z (r)}. 



(15) 



We see that the impurity-host coupling weakens rapidly 
for p }t 0.0. These figures show that, in 3D and for A = 2, 
the IBS is located at lujbs ~ 0.0, which is consistent 
with the results on ((M z ) 2 ) shown in Fig. 13. In Fig. 
13, for A = 1, we do not observe the development of 
a discontinuity for temperatures down to f3 = 32. For 
A = 2, we observe the development of a step centered at 
p 1=3 0.0, as (3 increases. For A = 4, a step discontinuity at 
p pa 0.3 is clearly observed. So, in 3D the IBS exists only 
for sufficiently large values of the hybridization matrix 
element V. This is consistent with the dependence of 
the IBS on the dimensionality in the {7 = case. These 




A 

s 

CN 



0.2 



0.1 



V 





i i i 1 

(3=16 




(c) A=4 




- - ^1=0.0 




" a--JX=-1.0 






a \ ■ a a - - 




* BJa-n-"'" 



10 

koR 



15 



20 



FIG. 11: (Color online) (Mi Mi) vs k R at @ = 16 and 
various p for hybridization (a) A = 1.0, (b) 2.0 and (c) 4.0 
for the two- impurity Haldane- Anderson model in the 3D case. 



results show that the dimensionality of the host material 
strongly influences the magnetic properties. 



The results presented in Sections 3 and 4 show that 
the density of states of the pure host at the gap edge and 
the value of the hybridization matrix element are crucial 
in determining presence and the location of the IBS. The 
IBS is in turn important in determining the magnetic 
properties of the systems when transition metal impu- 
rities are substituted into a semiconductor host. This 
means that the electronic state of the pure host material 
will also be crucial in determining the magnetic proper- 
ties. In the next section, we explore the consequences of 
a more realistic band structure for a GaAs host using the 
tight-binding approximation. 
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FIG. 12: (Color online) s(r) vs fcor at various j3 for A = 2 and 
for (a) /i = —0.1, (b) 0.0 and (c) 0.1 for the single- impurity 
Haldane- Anderson model in a 3D host. 



V. QMC RESULTS FOR THE TIGHT-BINDING 
MODEL OF A Mn IMPURITY IN GaAs 



In this section, we present QMC data obtained within 
the tight-binding model of a Mn impurity in GaAs, which 
was introduced in Section. 2.B. Here, we have per- 
formed the QMC calculations keeping all three of the 
Mn ti g orbitals in Eq. (1), hence this approach includes 
the multi-orbital effects except for the Hund coupling. 
In addition, we use a more realistic description of the 
semiconductor bands £k a for GaAs compared to that 
of Sections 3 and 4. Furthermore, the hybridization 
Vka.,g is determined by the tight-binding approach us- 
ing parameters consistent with photoemission measure- 
ments on Mn in GaAs, instead of being a free parameter. 
The following QMC results are obtained for the Slater- 
Koster parameters (pda) = —1.14 eV, (sda) = eV, and 
(pdn) = {pda) I '(-2.16). 

Figure 14 shows ((M z ) 2 ) and T\ versus p for the xy 
orbital near the top of the valence band for the single- 
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FIG. 13: (Color online) Impurity magnetic-moment square 
{(M z ) 2 } vs p, at various j3 for (a) A = 1, 2 and 4 for the 
single-impurity Haldane- Anderson model in a 3D host. 
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FIG. 14: (Color online) (a) <(A/ Z ) 2 ) and (b) T X versus p 
for a Mn xy orbital obtained using the tight-binding model. 
Here, the top of the valence band is located at /i = 0. 



impurity Haldane-Anderson model. Here, the develop- 
ment of a step discontinuity in the semiconductor gap is 
clearly seen as T decreases down to 180 K. At low T and 
for {pda) = —1.14 eV, the inflection point of the disconti- 
nuity occurs at 100 meV, which implies that ujibs = 100 
meV. This is in good agreement with the experimental 
value of 110 meV, especially if we note that the estimate 
of {pda) from the photoemission experiments is —1.1 eV. 

In order to understand how sensitively ujibs depends 
on the Slater-Koster parameters, we repeated these cal- 
culations for different values of (pda) keeping {sda) = 
and (pdn) = (pda)/ (-2.16). For (pda) = -1.4 eV, we 
found luibs ~ 300 meV, while for (pda) = —1.0 we ob- 
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FIG. 15: (Color online) Inter-impurity magnetic correlation 
function (Mf M|) vs it/a at various T for (a) /x = 0.0 eV, 
(b) 0.05 eV and (c) 0.10 eV for the two-impurity Haldane- 
Anderson model within the tight-binding model. Here, Aff is 
the magnetization operator at the Mn xy orbital. 



tained a more smeared discontinuity in ((M z ) 2 ) versus [i 
centered at /i w 0. 

Next in Fig. 15, we display (MfM|) versus the im- 
purity separation in lattice units, |R|/a, for the xy or- 
bital in the case of two Mn impurities substituted into 
Ga sites. Since we keep three orbitals at each Mn site, 
these calculations are more costly in terms of computer 
time. In this case, we present data for T = 1500 K and 
730 K. The Hirsch-Fye algorithm for the single and two- 
impurity Anderson model does not have the fermion sign 
problem. Hence, by using sufficient amount of computer 
time it is possible, in principle, to extend the calcula- 
tion of (MfM|) to lower temperatures. Even though 
the data on (Mf M|) are at high temperatures, the ef- 
fect of the IBS is observable. For T = 730 K and 
fi = 0.0, the FM correlations extend up to 1.5a, how- 
ever the FM correlations weaken as /i approaches loibs- 
We note that the spatial extent of the Mn-Mn interac- 
tion in low-density (Ga,Mn)As was also considered by 
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FIG. 16: (Color online) Impurity-host magnetic correlation 
function (M z m z (r)) vs r/a at various T for (a) /x = 0.0 eV, 
(b) 0.05 eV and (c) 0.10 eV for the single-impurity Haldane- 
Anderson model within the tight-binding model. Here, M z is 
the magnetization operator at the Mn xy orbital. 



using the tight-binding approximation and perturbative 
techniques in Ref. (29J. 

Finally, we discuss the correlations between the host 
electronic spins and the magnetic moment of the xy or- 
bital at the Mn site for the single-impurity case. In 
Fig. 16, we plot (M z m z (r)} versus the impurity distance 
|r|/a for various values of /i and T. Here, we denote the 
first-neighbor As site by As 1 , the second-neighbor Ga site 
by Ga 2 , and so on. These results show that the devel- 
opment of the AFM correlations between the magnetic 
moment of the xy orbital and the neighboring host elec- 
trons as T is lowered down to 360 K. We also observe 
the weakening of the AFM correlations as /i approaches 
^IBS- 
Using the tight-binding results for £k a and Vk Q ^ as 
input, in this section, we have performed the QMC sim- 
ulations to study the magnetic properties (Ga,Mn)As in 
the low-density limit. These results are similar to those 
presented in Sections 3 and 4 for the single-orbital case 
with simple band structure. In particular, we saw that, 
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by using a realistic model for the host band structure and 
host-impurity hybridization of (Ga,Mn)As, it is possible 
to obtain an accurate value for luibs- Such quantitative 
agreement supports the physical picture described in this 
paper for the origin of the FM correlations in DMS in the 
low-density limit. 

VI. DISCUSSION 

The computational approaches taken in this paper can 
be extended to study the case of finite density of Mn im- 
purities in GaAs. Here, we have reported results for the 
single- and two-impurity cases. However, the Haldanc- 
Andcrson model can also be studied for a finite den- 
sity of impurities using the Hirsch-Fye QMC algorithm. 
This type of investigation could shed light on the metal- 
insulator transition encountered in the DMS materials at 
finite density of impurities. In addition, it could help to 
understand the nature of the metallic state observed in 
(Ga,Mn)As for more than 2% doping. 

In the QMC results for the single-impurity case, we 
have observed the extended nature of the induced charge 
density and the spin polarization around the impurity 
site with a range lo determined by model parameters 
such as the hybridization, the bare band structure of the 
host, the Coulomb repulsion, etc. In the two-impurity 
case, when the polarization clouds overlap, the QMC re- 
sults showed that FM correlations develop between the 
impurities. Then, the interesting question is what hap- 
pens for finite density of impurities as the average sepa- 
ration of the impurities becomes comparable to l§. Does 
the Haldane- Anderson model capture the insulator-metal 
transition observed in (Ga,Mn)As? In addition, is it pos- 
sible to describe the electronic properties of the metal- 
lic state bordering the insulator-metal transition as that 
of a disordered valence band with weak correlation ef- 
fects? Or, does the metallic state retain more properties 
of the single-impurity case so that the IBS and the split- 
off states continue to play roles even though they emerge 
with the valence band? In this case, what is the nature 
of the magnetic correlations among the impurities? It 
would be interesting to study this problem even in the 
simplest case where each impurity has one orbital and 
the density of states of the bare host and the hybridiza- 
tion are constants as discussed in Section 3. By extend- 
ing the QMC calculations to the finite-density case, we 
think that it would be possible to address the question of 
whether the metallic state bordering the insulator-metal 



transition has unusual electronic properties. 



VII. SUMMARY 

In this paper, we studied the single- and two-impurity 
Haldane-Anderson models using the Hirsch-Fye QMC 
technique. Our purpose was to develop a microscopic 
understanding of the low-doping insulating phase of the 
DMS material (Ga,Mn)As. We first discussed the sim- 
ple case where each impurity consists of one localized 
orbital and the host band structure is described by two 
and three-dimensional quadratic bands. We have pre- 
sented QMC results on the inter-impurity and impurity- 
host magnetic correlations and the charge density around 
the impurity site. We observed that the presence and 
the occupation of the IBS is important in determining 
the magnetic correlations. We saw that long-range FM 
correlations, which are induced by the AFM coupling of 
the valence electrons to the impurity moment, develop 
between the impurities when the chemical potential is 
between the valence band and the IBS. We have also 
showed that, in the low-doping limit, the magnetic cor- 
relations between the impurities do not exhibit RKKY- 
type oscillations in a semiconductor. These calculations 
also displayed how the model parameters determine the 
magnetic correlations. 

In order to develop a more realistic model of the 
low-doping insulating phase of (Ga,Mn)As, we used the 
tight-binding approximation to map the system to the 
Haldane-Anderson Hamiltonian. In this approach, we 
used tight-binding parameters determined from the pho- 
toemission experiments and kept all three of the Mn ti g 
orbitals in the QMC calculations. The resulting value 
for loibs is close to the experimental value of 110 meV. 
We also observed that the magnetic correlations weaken 
as the IBS becomes occupied. These results are useful 
for developing a microscopic understanding of the low- 
doping insulating phase of (Ga,Mn)As. 
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